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Cross Reference To Related Application 

5 Reference is made to and priority claimed from U.S. 

provisional application Ser. No. 60/423,851, filed Nov. 5, 
2002, entitled METHOD AND APPARATUS FOR ESTIMATING THE 
IMPULSE RESPONSE OF A RADIO CHANNEL BASED ON A CALCULATED 
CORRELATION SEQUENCE. 

10 Field of the Invention 

The invention pertains to receiving a radio signal via a 
distorting radio channel, and more particularly to estimating 
the impulse response of the distorting radio channel so as to 
be able to adjust a receiver for receiving the radio signals. 

15 Background of the Invention 

As is well known, GSM (Global System for Mobile 
Communications) uses TDMA (time-division multiple access) in. 
which a radio receiver, such as the radio receiver 10 
indicated in Fig. lA, receives symbols communicated in one or 

20 more of the eight time slots of an 8-slot TDMA frame, the 

symbols being communicated as a radio signal by a radio 
transmitter (not shown) over a distorting radio channel (also 
not shown) . Distortion occurs for example because of 
multipath propagation which provides that at any particular 

25 time, the radio signal arriving at the receiver potentially 

includes contributions from the signal transmitted at several 
different earlier times. In general for a linear radio 
channel (or, approximately, for any radio channel) , if the 
transmitted radio signal is indicated as s(t) and the received 

30 signal as r(t) , then. 
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r(t) = c(t) ® s(t) + n(t) , ( 1 ) 

where c{t) is the impulse response of the radio channel, n(t) 

is the noise in the radio channel, and the operation ® 
indicates convolution . 

5 To account for distortion by the radio channel (or in 

other words, to account for an impulse response that provides 
other than an identity mapping of the transmitted to the 
received symbol sequence) , a receiver typically includes an 
equalizer, which is a module that in effect provides 
10 (approximately) a function h(t) that is the inverse (relative 

to convolution) of the channel impulse response c(t)) , so that 
the transmitted signal s(t) can be recovered from the received 
signal r(t) , i.e. the output s(t) of an equalizer is, 

Sit) = h(t)®r(t)^s(t) . (2) 

15 An equalizer is only one way of accounting for distortion 

caused by the radio channel . Another way is to use a module 
that implements the so-called Viterbi algorithm (well known 
in the art) , using the (estimated) channel impulse response 
and the received symbol sequence to arrive at (a best guess 

2 0 as to) the transmitted symbol sequence. The different ways of 

accounting for (i.e. in effect inverting) the impulse 
response c(t) of the radio channel are not the subject of the 
invention; the invention deals only with estimating c{t) . 

According to the invention and as in the prior art, the 

2 5 channel impulse response c(t) is modeled as a transversal 

filter, i.e. a tapped delay line, with each tap having a 
corresponding gain, generally a complex number; this follow 
from the widely used representation of signal as sequence of 
samples , 

3 0 =Z^f^A-i ' 
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where r^,Cp5^_. represent sampled versions of the signals 
r(t%c(t) and s(t) respectively. In the following, for 
simplicity, we assume operation at one sample per symbol, 
with an extension to a different sampling rate 
5 straightforward for those skilled in the art. A three-tap 

transversal filter is shown in Fig. IB. 

The aim of the invention is to estimate the gains c± of a 
transversal filter that represents the radio channel impulse 
response c(t) ; the number of taps of the transversal filter 
10 assumed to represent the radio channel impulse response c(t) 

is predetermined, based on the average duration of the 
channel impulse response and a judicious choice between 
accuracy and computational burden. 

Because the radio channel impulse response changes over 

15 time (since multipath, among other phenomena, is not 

constant, especially when either the radio receiver or the 
radio transmitter participating in a radio communication are 
mobile) , it is necessary to continually re-estimate the radio 
channel impulse response, i.e. to continually recalculate the 

20 gains Ci of a transversal filter assumed to represent the 

radio channel impulse response c(t) . To make such a 
recalculation possible, in the case of GSM, a burst (the 
transmission in a time slot) , which is 148 symbols long (each 
symbol representing one bit) , includes a training sequence X 

25 of twenty-six symbols, i.e. X = [x(0) x(l) x(2) ... x(25)] that 

in turn includes a sixteen- symbol correlation sequence Xc, 
i.e. Xc=[x(5) x(6) x(7) ... x (20)], sandwiched between two 
additional parts Xai Xa2 each five symbols long, as indicated 
in Fig. 2 and Fig. 3. (Xai = [x(0) x(l) ... x(4)] and Xa2 = 

30 [x(21) x(22) ... x(25)].) As shown in Fig. 4, the five symbols 

of the first additional part Xai ^^e the same as the last five 
symbols of the correlation sequence Xc, i.e. x(4)=x(20), 
x(3)=x(19) ... x(0)=x(16). Similarly, the five symbols of the 
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second additional part Xa2 are the same as the first five 
symbols of the correlation sequence Xc. (In the figures, the 
symbols of a burst are indicated as the corresponding bits.) 
Generally, the length of the additional part is chosen 
5 according to the number of taps of the transversal filter 

assumed to represent the radio channel impulse response c(t) , 
as will be clear from the description that follows. 

Training sequences are used not only in GSM, but also in 
other digital telecommunications systems and in essentially 

10 the same way as described above, i.e. by correlating a local 

replica with the received training sequence for various 
relative positions. (In such other systems, the number of 
bits per symbols is sometimes different than in GSM because 
different modulation schemes are used, the training sequence- 

15 length might be different than twenty-six symbols, the 

correlation length might be different than sixteen symbols 
and the memory of the channel impulse response (channel 
length - 1) to be estimated and consequently the length of 
the additional part might be different than five symbols.) 

20 The training sequence transmitted to the receiver 

detector module 11 (Fig. 1) is typically one of several 
different training sequences and the receiver detector module 
11 is provided with a code indicating the training sequence X 
being used by the transmitter. We indicate using Y the 

25 portion of the received sequence corresponding to the 

transmitted training sequence and refer to it as the received 
training sequence, writing Y=[y(0), y(l), y(2)... y(25)], where 

y(k) = ^CiX(k-i) , for 0<k<25 . 

i 

Similarly we refer to the portion of the received training 
30 sequence corresponding to the transmitted correlation 

sequence Xc as the received correlation sequence Yc. 
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To calculate the gains Ci, the receiver detector module 
11 performs sets of correlations of the received training 
sequence Y with a replica training sequence Xq (i.e. a copy of 
the training sequence known to have been transmitted) . The 
5 detector performs (calculates) correlations of the received 

training sequence Y with the replica Xo for typically six 
different relative positions of the replica Xq and the 
received training sequence Y, as described next; each 
correlation result yields a value of a corresponding gain Ci- 
10 For the GSM system, all the gains Ci can be estimated 

independently from each other because the used correlation 
sequences Xc included in the transmitted training sequence X 
(See Fig. 3) have ideal cyclic autocorrelation over a certain 
correlation lag. 

15 Fig. 5 indicates the correlations performed in a 

conventional channel impulse estimation (i.e. in a 
conventional calculation of the ci) . To perform the required 
correlations, a GSM receiver in essence slides the replica 
correlation sequence Xo,c over the received training signal Y 

20 one symbol at a time, starting with the received correlation 

sequence Yc and the replica correlation sequence Xo,c aligned. 
Thus, a correlation is performed for each of six different 
positions, each differing by one more symbol relative to the 
starting position. Each correlation provides the value for a 

25 corresponding gain ci (to within an overall normalization 

factor) of the transversal filter assumed to represent the 
channel impulse response c{t) , the correlation for the 
starting position yielding Co (see Fig. IB) , the next 
correlation for the next position yielding Ci, and so on. 

30 Thus, with a five-symbol (first) additional part, a six-tap 

transversal filter is used to represent the channel impulse 
response c(t) . 
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To obtain a more accurate determination of the channel 
impulse response c(t) , published European Patent application 
EP 0 701 334 A2, with Nokia Mobile Phones Ltd. as applicant, 
teaches a method of so-called double correlation, which makes 
5 use of the trailing additional part Xa2- As indicated in Fig. 

6, the double correlation method requires making two sets of 
correlations. A first set of six correlations is made in 
which, as in the conventional method, the received 
correlation sequence Yc is compared with a sequence consisting 

10 of the replica correlation sequence Xo,c and the replica first 

additional part Xo,ai/ starting with the received correlation 
sequence Yc aligned with the replica correlation sequence Xo,c 
and then shifting the received correlation sequence Yc one 
symbol at a time and repeating the correlation, until a total 

15 of six correlations are performed- Next, a second set of six 

correlations is made in which a different received 
-correlation sequence Yc' is used, consisting of only a portion 
Yc2 of the received correlation sequence Yc, namely the eleven 
symbols closest to the second additional part Ya2/ followed by 

2 0 the five symbols of the received second additional part Ya2/ 

as indicated in Fig. 7. It is useful to observe here, in 
order to understand the idea disclosed in the present 
invention, that in the absence of noise and under the 
assumption that the channel impulse response is effectively 

25 less than six symbol long, the symbols of the different 

received correlation sequence Yc' are the same as the symbols 
of the received correlation sequence Yc, although in a 
somewhat different order as indicated in Fig. 7. in the 
second set of correlations, the different received 

30 correlation sequence Yc' is compared with a sequence 

consisting of the replica correlation sequence Xo,c and the 
replica second additional part Xo,a2/ starting with the 
different correlation sequence Yc' aligned with the replica 
second additional part Xo,a2 and then shifting it one symbol at 
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a time and repeating the correlation, until a total of six 
correlations are performed. The results of the two sets of 
correlations (i.e. the gains ci) are then averaged (i.e. each 
correlation from the first set is averaged with the 
5 corresponding correlation in the second set) . 

The double correlation method thus requires twice as 
many correlations as the conventional method, and each 
correlation requires a relatively large number of arithmetic 
and logical operations by a processor. It would be 
10 advantageous to be able to somehow arrive at the same 

improved result achieved by the double correlation method but 
with substantially fewer processor operations than are 
required by the double correlation method. 



Summary Of The Invention 

15 Accordingly, in a first aspect of the invention, a 

method is provided by which a radio receiver, in receiving a 
signal transmitted over a radio channel, estimates the 
impulse response of the radio channel based on a received 
training sequence included in the received signal , the method 

20 including a step of performing a plurality of correlations of 

a correlation sequence derived from the received training 
sequence with a replica of the transmitted correlation 
sequence, characterized by a stepl) of calculating the 
received correlation sequence derived from the received 

25 signal based on averaging symbols of the received training 

sequence . 

In a second aspect of the invention, a method is 
provided for estimating the impulse response of a radio 
channel by which a radio receiver receives a received signal 
3 0 including a received training sequence for which the radio 

receiver knows a corresponding replica training sequence, the 
received training sequence including a correlation sequence 
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having a first end and a second and, and also including an 
additional part at the second end, with the additional part 
the same as a corresponding portion of the correlation 
sequence at the first end of the received correlation 
5 sequence, and likewise for the replica training sequence so 

that it includes a replica correlation sequence, the method 
characterized by: a sequence-calculating step, responsive to 
the received training sequence, of forming a calculated 
correlation sequence by averaging a predetermined number of 

10 symbols from the first end of the received correlation 

sequence with a predetermined number of corresponding symbols 
from the additional part at the second end of the received 
training sequence; and a correlating step, responsive to the 
calculated correlation sequence, of performing a set of 

15 correlations of the calculated correlation sequence with the 

replica training sequence, the set of correlations including 
a first correlation in which the calculated correlation 
sequence is aligned with the replica correlation sequence and 
including subsequent correlations performed with the 

2 0 calculated correlation sequence shifted for each next 

correlation by one or more symbols from the position in the 
immediately preceding correlation, so as to provide 
information useful in estimating the channel impulse 
response . 

25 In a third aspect of the invention, a receiver is 

provided, characterized in that it is operative according to 
either the first or second aspect of the invention. 

In a fourth aspect of the invention, a telecommunication 
system is provided, including a base transceiver station and 

3 0 a user equipment, both of which include a receiver, 

characterized in that both receivers are operative according 
to either the first or second aspect of the invention. 

In a fifth aspect of the invention, a computer program 
product is provided comprising: a computer readable storage 
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structure embodying computer program code thereon for 
execution by a computer processor in a receiver, with said 
computer program code characterized in that it includes 
instructions for performing steps of a method according to 
5 either the first or second aspect of the invention. 



Brief Description of the Drawings 

The above and other objects, features and advantages of 
the invention will become apparent from a consideration of 
the subsequent detailed description presented in connection 
10 with accompanying drawings, in which: 

Fig. lA is a block diagram showing some components of a 
receiver of a cellular telecommunications device, including 
a detector implemented to perform the method provided by 
the invention. 

15 Fig. IB is a block diagram of a transversal filter (a 

tapped delay line) assumed to represent the radio channel 
providing the radio receiver with a received radio signal 
differing from the transmitted radio signal because of 
distortion by the radio channel (modeled by having non-zero 

2 0 gains for other than only the first tap of the transversal 

filter) . 

Fig. 2 is a schematic illustrating the frame structure 
of a frame used in time division multiple access (TDMA) by 
some cellular telecommunication systems (such as GSM) , and 
25 also showing the composition of a so-called normal burst 

(the transmission in one time slot of a TDMA frame) , 
including a training sequence. 

Fig. 3 is a schematic illustrating the composition of a 
training sequence . 

30 Fig. 4 is a table showing the bits (corresponding to 

symbols of the training sequence via a mapping such as a 
mapping in which bit value 1 corresponds to symbol value 1 
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or some other symbol value, and bit value 0 corresponds to 
symbol value -1 or some other symbol value different from 
the value corresponding to bit value 1) in each of the 
eight different training sequences used in GSM. 

5 Fig. 5 is a schematic illustrating the correlations 

performed by the conventional method for using a training 
sequence (in GSM and some other telecommunication systems. 

Fig. 6 is a schematic illustrating the correlations 
performed by the so-called double correlation method of the 
10 prior art. 

Fig. 7 is a schematic illustrating the composition of 
the two different correlation sequences Yc and Yc' used in 
the so-called double correlation method of the prior art. 

Fig. 8 is a schematic illustrating the composition of 
15 the correlation sequence Yc" used in the invention, which 

is a calculated correlation sequence. 

Fig. 9 is a schematic illustrating the correlations 
performed according to the invention, using the calculated 
correlation sequence Yc" - 

2 0 Fig. 10 is a flow chart indicating the method of the 

invention. 

Best Mode For Carrying Out The Invention 

The invention is described below making reference in 
particular to GSM. It should be understood however that the 

25 invention can be used in any telecommunication system using 

TDMA and using training sequences constructed similarly to 
those used in GSM, i.e. so as to have a correlation sequence 
and two additional parts, one on either side of the 
correlation sequence, with the correlation sequence including 

30 fragments that are the same as the two additional parts . 
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Referring now to Fig. 1, a receiver of a 
telecommunications system 10 is shown according to the 
invention as including a detector 11 for demodulating a 
received signal of a stream of symbols modulating a carrier 
5 wave, a module for estimating the channel impulse response 

c{t) (i.e. for determining gains ci of a transversal filter 
representing the channel impulse response c(t)) , an equalizer 
or other means of accounting for the channel impulse response 
c{t) (such a module implementing the Viterbi algorithm) and so 

10 determining a transmitted symbol sequence corresponding to a 

received symbol sequence, given the estimate of the channel 
impulse response c{t) , and a symbol mapper for mapping the 
received symbols to corresponding bits. The received symbol 
stream is a so-called normal burst (the content of a time 

15 slot of a TDMA frame) . In at least some of the bursts, a 

(transmitted) training sequence X is conveyed and received as 
a received training sequence Y, providing the receiver with a 
means for estimating the channel impulse response c{t) in 
terms for example of gains ci of a transversal filter, as 

20 discussed above, by performing correlations of the received 

training sequence Y with a replica Xq of the transmitted 
training sequence X (the replica being stored in the detector 
11) . Since the transmitted training sequence X is only one of 
various different training sequences available for use, the 

25 detector is either provided with a code indicating the 

particular sequence being used, as shown in Fig. 1, or else 
determines the training sequence being used by examining the 
received bit stream. Ultimately, the detector provides a bit 
stream (only the data bits of the received bit stream) to 

30 (typically) a decryption module 13, which, in combination 

with other modules (not shown) , further processes the 
detected bits, and then a channel decoding module 14 decodes 
the further processed bit stream (removing the redundancy 
added to allow error detection and possibly error correction) 
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and a speech decoder 15 and a digital to analog converter 16 
process the output of the channel decoding module 14 so as to 
provide a speech signal to a microphone 17. 

The above description refers to voice transmission, but 
5 it is to be understood that the invention is also of use in 

case of data transmission. 

Referring now to Fig. 8 and also to Figs. 3 and 7, 
instead of using either a portion of the received correlation 
sequence Yc as in the conventional method, and instead of 

10 using a different correlation sequence Yc' (=[Yc,2 Ya2 ]/ as 

shown in Fig, 7) used in the double correlation method of the 
prior art, the invention calculates a received correlation 
sequence- -called here the calculated correlation sequence and 
indicated by Yc"--so as to include the eleven- symbol portion 

15 Yc2 (consisting of the eleven symbols closest to the second 

additional part Ya2) of the received correlation sequence 
concatenated with a five -symbol string Yd' that is an average 
of the five symbols Yd (the five symbols closest to the first 
additional part Yai of the received training sequence Y) and 

20 the five symbols of the second additional part Ya2 of the 

received training sequence Y. Thus, the calculated 
correlation sequence Yc" is given by, 

Yc" = [(Yci+Ya2)/2 Yc2 ] = [Yci' Yc2 ] , 

where the square brackets [] indicate concatenation of the 
25 sequences inside the brackets, i.e. Yd' and Yc2- 

As a result of the averaging process used for the first 
portion of the calculated correlation sequence Yc" , the noise 
in the calculated correlation sequence is diminished. It will 
also be evident for those skilled in the art that there are 
30 alternative ways to obtain the calculated correlation 

sequence (correlated by the proper portion of reference 
training sequence) which would provide the same estimation of 
the gains ci, as for instance. 
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alternative — [Yc2 (Yci+Ya2)/2] = [Yc2Yci']. 

Referring now to Fig. 9, instead of performing two 
correlations as in the double correlation method of the prior 
art, the invention performs a single correlation as in the 
5 conventional (prior art) method, but uses the calculated 

correlation sequence Yc" instead of the received correlation 
sequence Yc- The result is the same as what is provided by the 
double correlation method but requires approximately one half 
the number of processor operations (i.e. approximately the 
10 number required to perform one set of six correlations) since 

the number of processor operations required to provide the 
calculated correlation sequence Yc" is substantially less than 
what would be required in performing a complete additional 
set of six correlations. 

15 Referring now to Fig. 10, the invention is shown as 

providing a method including a first step 101 of forming the 
calculated correlation sequence Yc" as indicated above, or 
more generally, by averaging a predetermined number of 
symbols from the second additional part Ya2 of the received 

2 0 training sequence Y with a corresponding predetermined number 

of symbols from the end of the received correlation sequence 
Yc away from the second additional part Ya2 ■ Then in a next 
step 102, the correlations (typically six) are performed as 
in the conventional method but using the calculated 

25 correlation sequence Yc" in place of the received correlation 

sequence Yc, i.e. correlations of the calculated correlation 
sequence Yc" with the replica training sequence Xq starting 
with the calculated correlation sequence Yc" aligned with the 
replica correlation sequence Xo,c and beginning each 

30 subsequent correlation shifted one symbol toward the first 

additional part Xo,ai of the replica training sequence Xq . 

It should be understood that the above described process 
is also of use in case of a modulation system in which more 
than one bit is conveyed per symbol. 
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The invention has been described in terms (primarily) of 
the steps of a method. The invention also comprehends an 
apparatus- -e . g . the detector 11 (Fig. l)--for performing the 
above described steps. For each step described above, there 
5 can be a corresponding module of the apparatus, although it 

is also possible for the functionality for performing more 
than one of the above-described steps to be incorporated into 
a single module. Such modules may be implemented as hardware, 
or may be implemented as software or firmware for execution 

10 by a processor. In particular, in the case of firmware or 

software, the invention is provided as a computer program 
product including a computer readable storage structure 
embodying computer program code- -i.e. the software or 
firmware- -thereon for execution by a computer processor 

15 provided with e.g. the receiver 10 (Fig. 1) of which the 

apparatus such as the detector 11 (Fig. 1) may be a 
component . 

It is to be understood that the above -described 
arrangements are only illustrative of the application of the 

20 principles of the present invention. In particular, the 

receiver 10 (Fig. 1) can be the receiver in a mobile station 
having a radio channel between its antenna and a base station 
antenna, or it be the receiver in the base station (as 
already mentioned) , or it can be a receiver in a telephone 

2 5 connected via a cable to a remote transceiver with a radio 

channel between the transceiver antenna and the base station 
antenna. In addition, the training sequence need not 
necessarily be in the middle of a burst. Thus e.g. in the 
digital radiotelephone system used in the United States, a 

30 TDMA frame consists of six time slots of 162 symbols each, 

and a burst transmitted in a time slot from the base station 
to the mobile station always begins with a synchronization 
burst of 14 symbols (typically representing 28 bits) , which 
is used as a training sequence. Also, the additional part 
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need not be on both sides of the correlation sequence, as in 
the GSM system, but instead there could be only one 
additional part which either before of after the correlation 
sequence- Finally, the symbol sequences in the training 
5 sequences need not necessarily be the same in both directions 

(from base station to subscriber, and from subscriber to base 
station) . Numerous modifications and alternative arrangements 
may be devised by those skilled in the art without departing 
from the scope of the present invention, and the appended 
10 claims are intended to cover such modifications and 

arrangements . 
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